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Abstract 



The rate of 5-violation in the standard model at finite temperature is closely 
related to the diffusion rate F of Chern-Simons number. We compute this rate for 
mil ~ in the classical approximation in an effective SU(2)-Higgs model, using 
Krasnitz's algorithm. The parameters in the effective hamiltonian are determined 
by comparison with dimensional reduction. In the high temperature phase we 
find T /V{a\YTY ~ 1? neglecting a finite renormalization. In the low temperature 
phase near the transition we find the rate to be much larger than might be 
expected from previous analytic calculations based on the sphaleron. 
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1 Introduction 



Fermion number is not conserved in the standard model of elementary particle 
interactions, because of the anomaly in the divergence of the B + L current. This 
leads to a rate F of fermion number violation which plays an important role in 
theories of baryogenesis [||, H. In the approximation of taking into account only 
the SU(2) contribution to the anomaly, F can be identified from the diffusive 
behavior of the 'topological susceptibility' at large (real) times [^, 

rt = (g'(t))T, t^oo, (1) 

Q{t) = j^dx' j d'xqix), q = ^^F^^F^^\ (2) 

Here F is the SU(2) field strength and the bracket in (|l|) denotes the expectation 
value at temperature T. 

For temperatures below the electroweak phase transition, analytic calculations 
based on the sphaleron solution lead an expression which may be summarized in 
the form |, |, § 

— = K{awTY, (3) 




where V is the volume, aw = 5'^/47r is the SU(2) gauge coupling, X^/g^ ~ X/g^ ~ 
mjj/8m'^ is a coupling ratio in the 3D dimensionally reduced theory and Es{T) 
is the temperature dependent sphaleron energy. The effect of the fermions enters 
mainly in the relation of the 3D couplings to the original 4D standard model 
(see also [§]). 

For high temperatures above the phase transition, dimensional reasoning leads 
to the form (^ with k roughly temperature independent [|l], || (a logarithmic 
dependence may still be present). An analytic estimate k ^ 0.01 was carried out 
in §. 

The rate F defined by (P is a nonperturbative quantity and it is important 
to check analytic calculations by numerical simulation. This is well-known to be 
difficult since the effective Boltzmann factor in real time processes is complex: in 
the temporal gauge Aq = we can write 

Q{t) = iVcs(t)-iVcs(0), (5) 
{Q'it))r = Tre-^/ne^''*^cs(0)e-^^*-iVcs(0)]VTre-^/^, (6) 

where A^cs is the Chern-Simons number and the trace is over states in Hilbert 
space satisfying the Gauss constraint. To cope with the complex weights a 
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classical approximation has been introduced in |]T0| and tested on the abelian 



Higgs model in 1+1 dimensions |Tl[], using the microcanonical ensemble. Subse- 
quent computations used Langevin methods |T2|, O and the canonical ensemble 
[ffl, IT^, |TB|. In these computations the quantum mechanical expectation value 



is replaced by a classical expression, 



where ip and tt denote generic canonical variables and S{G) enforces the clas- 
sical Gauss constraint. Furthermore, ip{t) and 7[(t) are solutions of Hamilton's 
equations with hamiltonian i^eS and initial conditions if{0) = (f, 7r(0) = vr. The 
classical partition function is still a functional integral over all v'(x) and 7r(x) and 
needs regularization. In the numerical simulations this is provided by a lattice. 

In the abelian-Higgs model the numerical results agree with the analytical 
sphaleron-type calculations in the low temperature regime [|14|, |15|. In the high 



temperature regime the rate appears to follow the 1+1 dimensional analogue 
of @ with temperature independent k, which depends, however, on the lattice 



spacmg 



(see also ||T5|)- As emphasized in [T^, such a lattice spacing de- 
pendence should be cancelled by a corresponding dependence in the hamiltonian 
Hes, which is to be an effective hamiltonian appropriate for the lattice used. 

In the realistic 3+1 dimensional case of the SU(2)-Higgs model, exploratory 
simulations have been carried out in [|T^. The incorporation of the Gauss con- 
straint has been a source of uncertainty in the nonabelian case. Recently, a 
solution to the problem has been presented |18| and first results have been ob- 
tained for K in the pure SU(2) case |]2D|. Here we study the full SU(2)-Higgs 
case using the algorithm and implementation of ref. |]T8[. As this work was being 
completed a paper appeared describing simulations with a chemical potential for 
the Chern-Simons number, using a different implementation of the Gauss law 



IS 



In sect. 2 we discuss the classical approximation in general terms and em- 
phasize the relation to dimensional reduction (DR). In sect. 3 we describe the 
determination of the parameters in the effective hamiltonian for the SU(2)-Higgs 
model by comparison with the dimensional reduction results of refs. 0, |2^. The 
connection with dimensional reduction is also discussed in refs. [|], In 
sect. 4 we go through some conventions used in the numerical implementation 
and in sect. 5 we present results for the rate F in a temperature region around 
the phase transition. Sect. 6 contains our conclusion. Some details are delegated 
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to the appendices. 



2 Classical approximation and dimensional re- 
duction 

A classical approximation is expected to be valid if expectation values are domi- 
nated by states with high mode occupation numbers. This will occur for tempera- 
tures higher than any mode energy, which suggests introducing a cutoff A to keep 
these energies bounded from above. We can imagine integrating out the modes 
on spatial momentum scales larger than A, leading to an effective action S'efr for 
the remaining variables. If this effective system is weakly coupled, such that a 
description in terms of 'modes' applies, with mode energies < a/ (mass)^ + A^, 
we may expect classical behavior for A ^ T. 

For gauge theories a separation in low and high energy modes is notoriously 
difficult. However, sensible approaches do exist, for example real space renormal- 
ization group transformations in lattice formulations , or methods as outlined 



in 



For definiteness we assume A to be represented by a spatial cubic lattice 
with spacing a = vr/A. Keeping time continuous seems artificial and perhaps 
we should also coarse-grain in time. For slow processes a derivative expansion 
may be adequate and keeping only two time derivatives we can go over to a 
canonical formalism with an effective hamiltonian H(,g{(f,TT). The procedure is 
to be applied to observables O which are dominated by low energy modes which 
can be expressed in the effective variables in a practical way, O Ocs{'^, tt). 

Under weak (effective) coupling conditions this hamiltonian may then 
be approximated by a classical form in which the parameters depend on A and 
T. To avoid strong coupling in a nonabelian gauge theory, the scale A should 
presumably not be too low. The procedure thus leads first to an expression of 
the form (^), but with H —>■ ifgff on a spatial lattice. This quantum mechanical 
expression is then further approximated by the classical expression (|^. 

In non-gauge theories the hamiltonian Hes has the typical form 

i7eff = $^^vr(xf + (9) 

X 

where W includes the spatial gradients. Performing the gaussian integration 
over the canonical momenta vr in the partition function gives a reduced partition 
function which has a dimensional reduction form 

Zdr = exp[-SMv>)], Sj,M = W{^)/T. (10) 
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Correspondingly, the static (equilibrium) expectation values in the classical effec- 
tive theory have the dimensional reduction form. So we may identify the classical 
approximation of the effective theory with a dimensionally reduced theory on a 
lattice. 

The dimensional reduction method can be formulated as a matching procedure 
[ P^ , in which the parameters in a local DR action are calculated in perturbation 
theory to give the same physical results as the original 4D theory, approximately. 
The error in the approximation is partially due to the neglect of nonlocal terms 
which appear first at the two loop level p6| . 

Now the cutoff A in the DR theory is interpreted as an ultraviolet cutoff, which 
is to be removed to infinity while tuning the DR parameters in the appropriate 
way. However, it seems undesirable to remove the regulator in the reduced theory 
completely, since this would imply tuning to a second order phase transition, with 
accompanying loss of free parameters due to universality. For example, in the 
SU(2)-Higgs theory the ratio mj^/mvy would be fixed at the critical point. Such a 
problem does not enter in the approach based on integrating out high momentum 
modes, in which the cutoff is not be sent to infinity, which implies of course some 
nonlocality. 

But this method has its own awkward features when A ^ T. Then all modes 
up to A are highly excited and the lattice nature of the system will be ubiqui- 
tous. To obtain cutoff independent results under such conditions we may have 
to take into account the regularization dependence of the effective observables. 
But the calculation of the effective hamiltonian is already difficult and has to our 
knowledge not been done so far, let alone the obervables. 

However, the phenomenon of dimensional reduction suggests considerable 
freedom in the choice of A, including A ^ T, provided the observables are dom- 
inated by low momenta. Using such cutoffs A ^ T, regularization artefacts get 
suppressed and the latticized classical form of the observables will suffice (tak- 
ing into account necessary renormalizations) . Most importantly, we are lead to 
identify the potential energy part of the effective hamiltonian with S'dr- 

The initial conditions in (|^) now correspond to the DR action whose renor- 
malization properties are well understood. For short times one expects the real 
time correlators to have similar renormalization properties because Hamilton's 
equations imply a continuous dependence on the initial conditions. An example 
of a one loop perturbative calculation in scalar field theory can be found in ref. 



2^. For large times the dependence on initial conditions is governed by Lya- 
punov exponents, which should be renormalized quantities (i.e. of the order of 
a physical mass scale and not of order of the cutoff) in order to have renormal- 
izability of the real time correlators. Evidence for this is provided in ref. [13], 
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where Lyapunov exponents were measured in SU(2) gauge theory. The difference 
between two neighboring initial conditions a and h was characterized by electric 
and magnetic 'distances' De,m, 



De oc / d^x \El^) - E, 



Dm oc / (Tx \Bl) - B. 



(11) 



(in continuum notation). Assuming ultraviolet divergences to cancel in these dif- 
ferences, the Lyapunov exponents extracted from their exponential time behavior 
may be expected to be renormalized quantities. In ref. it was foundQ that 



De and Dm lead to the same Lyapunov exponent A = 0.30 ^f^T. 



3 The effective gauge-Higgs hamiltonian 

We now look more closely at the SU(2)-Higgs theory with classical action 

j d'x[^^F^,F^^'^ + {{d,-iAM\d^-iA^)^ + ^^^^^ + X{^'^^fl (12) 

where (p is Higgs doublet and = A'^Ta/2, with Tq, the three Pauli matrices, 
a = 1,2,3. Our basic assumption is that a reasonable approximation to the 
effective action is provided by a lattice version of (O): 



'S'ofr 

LeS 

W 



dx" Leff, 



(13) 



■ Tr {DoUm^yDoUmx + — {{do - iAo^)(f^y{do - iAo^)if^ 



^^Tl{l- Umnyi) + (I^mV^x)''' A 



doUr, 



^^Ox^^mx + ^f^mx^Ox+mj -^mV^x — Umx^y 



(14) 



(15) 
(16) 



where we used lattice units a = 1. The Umvi ^ exp(— iaAmx) are the parallel 
transporters at x for directions m = 1,2,3 and t/mnx is the transporter around 
a plaquette. The couplings ze, ^tt, qIs, /^es ^^"^ -^es depend on the temperature, 
lattice spacing and the renormalized couplings (7^, /x^, A in some renormalization 
scheme. We have kept the standard normalization for the spatial gradients, which 



•^In the notation introduced in (29),(p3|) the result in was a\ — 1.20//3, where a is the 
lattice distance and /3 = A/g-^aT. 
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defines g'^^ and the other couphngs. The kinetic terms have renormalization 
factors ze and z^^, which will be of the form ze.-w = 1 + 0{g'^, A). They influence 
the time scale relative to the momentum scale. The constant e adjusts the energy 
density and plays no dynamical role in our theory (it takes care of the Rayleigh- 
Einstein- Jeans divergence). 

The action is invariant under gauge transformations Q:x_[x^) acting in the usual 
way. A canonical description is straightforward in the temporal 'gauge' Ao = 0, 
leading to the hamiltonian 

^eff = E ( Y 9eEK.K. + ^.4^x) + W, (17) 

X 

where E^^ is the left translator E^^{L) or right translator E^^{R) (cf. appendix 
A for more details), and vr^ is conjugate to yjx- The residual invariance under 
static gauge transformations is dealt with by imposing the Gauss constraint. In 
the quantum theory G^lphys) = 0, with 

= ^ = i'D'^E-^+f^-), (18) 
D'mEm^ = y^}E^^{L) - E^^_^{R)], (19) 

m 

■Oct • f ''"a ■ i Ta /nn\ 

Jx = ^'/^LyTTx-Ky </'x- (20) 

In the classical approximation, the partition function entering in (|^ reads for the 
SU(2)-Higgs theory (cf. appendix A) 

Z = J DEDTTDUDip[Yl5{G'^)] expi-H^s/T). (21) 



If we now write 



= / ^ exp(.A^,GS), (22) 



carry out the integrations over E and vr and subsequently rescale A^^ — > ^/zeAq^/T, 
we arrive at the dimensional reduction form 



ZuR = J DAoDUDipexpi-SuR), (23) 



with 



^DR = ^ + ^E[^^-^Oxi^-^Ox + ^^e'ff¥'i¥'xAox^^x], (24) 

5'ofr-' X 

-Dm^Ox = -Ra/3(t^mx)^0x+m ~ ^Ox (25) 
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{Rap{U) is the adjoint representation of U). This fits into the general DR form 
used in 1^1], ^ in numerical simulations, 



5, 



DR 



/^G 5Z (1 - ^ Tr U^n.) - /5h 5^ i Tr ^lUm^^ 

xm<n xm 

+ 5^ i Tr + Tr <^l<^^f 

X X 

+ /5g 5^ ^ Tr (A^^ - Aoxf/L^ox+mf/, 

xm 

-/3,^5^lTrAL + /5f$:(^TrA: 



.2 ^2 
^OxJ ? 



+ ^yiTrA2,iTr<|.t$. 

X 

(26) 



with 



R ^eff 



Ze_ _ . 



f^2 = Pi 



PgPr 




2 2 



2(1 - 2/3fi - 3Ph) 
(3h 



(27) 
(28) 



(we have redisplayed the lattice distance a). It is clear that in customary DR 
notation g^^T = g^, XesT = A3, for lattice spacing a — in perturbation theory. 

The parameters Pg etc. are related to the 4D couplings, lattice spacing and 
temperature, depending on the 4D model [^. For the SU(2)-Higgs mo 

4 



g^aT 



m 



H 



4r2 



-A 



G5 



1 

~2 

9' 
+ — 
2 

5 4 



1 1 



P(3h 9 



47r/3c 



^-^lln 
9 27 



^ + ^ + ¥^-27^ 



9' 



16 ^ 16 ^ 167r2 V 96 



149 3p 

+ 35 



3 g^g. 



G 



lOS, A 



487r2 



2 ' 



(29) 



(30) 
(31) 
(32) 



with S = 0.252731, r] ^ 2.18, r/ ^ 1.01, ^ 0.44 IgT], |2g. The 4D couplings 
A, and MS* scheme quantities at the scale /i^ = 47rTe ^ 7T [pT[|. 

■^These formulas taken from H, E3] contain an error which is negligible for our purpose (of. 



the comments below eqs. (142), (143) in 23 
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The mass ratio m'j^/m'^ is approximately equal to its zero temperature value. 
We also have p ~ SX/g"^ ~ SXcs/Qes, with small corrections of order g'^. 

The term in the DR action is negligible in practice {g ~ 2/3). The 
coefficient ratio ^2/^0 is the Debye mass m|, = 5g'^T^/12 corrected with a coun- 
terterm -5^2TS/12. Setting (3^ to zero means that the Debye mass generated 
by ( pT| , piD will not have its standard perturbative value, but the accuracy of this 
perturbative value is unclear anyway. Our numerical results for the transition 
temperature to be described later indicate that the couplings in the effective 
action are not very important. 



The relation ze/zt^ = 1 is an approximation used in |22|; the corrections 



are small (cf. eq. (18) in ze/z^, = h^/ {g1/ 4)). With the approximation 
ze = = z, this parameter only affects the time scale. It takes care that the 
velocity of light c = 1. In the following we shall absorb z in the time scale (i.e. 
set z = 1), which means that c 7^ 1 in our units. However, we expect that 
z = 1 + 0(5'^, A) is close to 1 anyway. 



4 Numerical implementation 

In the classical approximation the rate F can be computed by numerical sim- 
ulation according to (|^). We use the algorithm and numerical implementation 
offered in ref. |T^]. For clarity we record the conventions used in [rH|, obtained 
by extracting an overall factor 4/5fgg from the effective lagrangian ([T5|) and going 
over again to canonical variables (cf. appendix A). Indicating the quantities of 
ref. |jl8| by a bar we have 



eflf 



T 
H 



Y,[\ E'^A. + + J] (1 - ^ Tr 



m<n 



with 



V . Aeff _2 2 2 9eS 

X = 4—, V =-a /ieff^^- 

9cS °^efi 



2\2l 



Hamilton's equations follow from the nontrivial Poisson brackets 

Our investigation has concentrated on the case mn ~ wiiy. Using eq. 
fix the parameter A = 1/2 since 2 A = 8Xcs/g^g ~ S-^/S'^ ~ ''^'h/t^w ~ 



(33) 

(34) 
(35) 

(36) 

]) we 
= 1. 
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The temperature and lattice spacing dependence of A can be neglected. The 
quantities Pr and v"^ are then given in terms of Ph and Pg by (|27|) and (|35|), in 
particular 



_2 2 ^ pPh 1 



Eq. ( pUD gives an analytic constraint on the lattice spacing and temperature 
dependence of v^. It is instructive to rewrite (^) in the form 

-2 2 2 V (3 + p 3 + p2S\ 2^2 V 1 /149 + 9p 

4Aw = a - - ^ — ] aT - [-g' 



X) V 18 3 aT J V2 / 87r2 ^ 486 

27 2 27 y ' ^ ^ 

which shows how depends on a and T if we neglect the temperature variation 
of niH- The latter runs rather slowly with T at /ij^ ^ 7T and T near Tg. We also 
neglect the running of and fix this coupling as (7 = 2/3; then aT = 9/(3. The 
temperature in units of tuh now follows from (^) or (^). 

We also follow in the lattice implementation of the topological charge 
density, 

k 

(39) 

with k,l,m = l,2,3(cycl) (the sum is over identically oriented Wilson loops in 
the plane perpendicular to direction k along the four plaquettes starting and 
ending at x). The normalization of is guided by the continuum limit a ^ 0. 
The magnetic field strength Fj/^^s) is identified from any of the four plaquettes in 
( |39D in the usual way and the electric field strength is identified from — 2£'^(x) — *• 
9oA^(x) = by Hamilton's equations. 

In the context of zero temperature QCD this lattice version of g(x) is some- 
times called a 'naive' implementation (see p9| for example for a review). As there 
is no other gauge invariant pseudoscalar field of dimension four to mix with, one 
expects that a finite multiplicative renormalization is sufficient to obtain the de- 
sired observable, 

g(x) = K^ql (40) 

This renormalization may be substantial. In perturbation theory Kq is largely 
due to leaf (or 'tadpole') diagrams, resulting from the compact nature of the 
lattice gauge fields, but the perturbative value is often unrehable. In the QCD 
context values for Kg of order 5 have been reported in ref. pO[ (where Kq is 
denoted by Z~^). Ideally one would like to use fermions to measure directly 
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N 




-vl 


Tc/niH 


12 


8 


0.3487 


0.27894 


1.91 


12 


20 


0.347733 


0.26295 


2.14 


12 


24 


0.34772 


0.26273 


2.15 


20 


32 


0.34173 


0.15597 


2.15 



Table 1: Critical values of (/9,f) based on the DR data for {Pg^ f^n) presented 
in ref. [g2]. 



how they are affected by changes in the topology of the gauge field, as in the 
'fermionic approach' to topological charge in QCD PT| , For an interesting 
study of chirality and zero modes in the present context see ref. [Bl]. We have 
not calculated k„ yet. 



5 Numerical results 

Numerical results for equilibrium quantities using dimensional reduction have 
been obtained in ref. ||22|. We use the same parameters at the phase transition as 
in this paper (except that /3^4 = 0, of course). Using the critical values {Pq, Pjj) 
at various lattice sizes from |^ gives v'^ according to eq. (|37|), shown in 



Table |. 

Since our effective hamiltonian implies vanishing couplings 4 '^^ ^^e Aq field 
in the associated DR theory, we want to see if this has an effect on the critical 



temperature. In ref. |^ Pg has been fixed while Ph "was varied to find the phase 
transition. Instead, we have kept v"^ fixed at v"^ and varied (3. An example is given 
in Fig. |I|, where v"^ corresponds to Pq = 12 in We see that Pc ~ 12 as well. 
From similar results for other observables at various volumes and pairs {Pq, Pjj) 
we have concluded that in all cases the critical Pc was approximately equal to the 
input Pq used for the calculation of (we have not made a precise determination 
of the transition values Pc, using e.g. finite size scaling). Hence, the fact that the 
effective hamiltonian implies /5^4 = does not seem to be important for the 
couplings used in our study, at least not for the value of Pc- In the following we 
set Pc = Pg, as indicated in Table |l|. 

Different values of Pc correspond to different lattice spacings a = 4/g'^Tp. As 
a consistency check we may compare the a-dependence of Vc obtained from the 
measured values of Pjj in ref. ||2^ with that predicted by eq. (^). The two pairs 
{P"^, N) = (12, 20) and (20,32) are suitable for such a check as they correspond to 
approximately the same physical lattice size L = Na: 32 ^ (20/12) * 20 = 33.3. 
Starting from the measured value of = 0.26295 at {Pc, N) = (12, 20), eq. (^) 
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5 10 15 20 



Figure 1: The correlation v^x^/mxV^x+m/lv'xl |v^x+m| as a function of P for N = 24 
at fixed = = 0.263. 

gives v'^ ^ 0.162 at {(3c,N) = (20,32) (using g = 2/3 and p = 1), which is close 
to the measured value ~ 0.156 in Table |I]. Neglecting finite size and additional 
lattice spacing effects, the critical (3c and v'^ values lead via (|38|) to the ratios of 
Tc/mn shown in the right column of Table 0. 

Our results for the Chern- Simons diffusion rate do not incorporate yet the 
renormalization factor Kg. An illustration of the effect this can have is in Fig. ^ 
which shows the Hamilton evolution of Qhii) = Jq dx^ Sx^x starting from a 
configuration equilibrated at /3 = 36, with v"^ = 0.1564 corresponding to (3c = 12, 
v'^ = 0.279 via (|3^). This is a rather large (3 value corresponding to a low 
temperature T ^ Tc/3 (neglecting finite size effects). We see fluctuations about 
classical vacua separated by AA^'cs ~ 0.75. With Hg = 1/0.75 = 1.33 we would 
get AiVcs ~ 1- Writing this as Kg ^ 1 + 12//? would give Hg ^ 2 at (3 = 12, which 
indicates that the Kg correction can be substantial indeed. 

The fact that the difference of Chern-Simons number between two adjacent 
vacua in Fig. |^ is near 1 also provides an indirect check on the normalization of 
the lattice form of q^. This is relevant in view of the fact that the numerical rate 
turns out to be very different from expectations based on sphaleron calculations. 

The volume in Fig. |^ was chosen small (in addition to large (3) in order to 
suppress fluctuations of different portions of the volume, which would obscure 
the classical vacua. Fig. |^ shows a typical example of the time evolution of Qhi^) 
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Figure 2: Hamilton evolution of = A/'cs(^) -^cs(O) 

on an 8 lattice starting 
from an equilibrium configuration at /? = 36, v'^ = 0.1564. 



for a 24^ system at /? = 13. Fig. ^ shows the corresponding average {Q\{t)) over 



33 configurations from the canonical ensemble. Following ref. |^0| we shall reduce 
the errors by taking also a microcanonical average, by first averaging over the 
Hamilton part of the trajectories and then over the intitial conditions obtained 
with the Langevin algorithm: 



average 





fto+t 


/ dto 


/ dt' 


h 


J to 



dt'J2^tit') 



<t<tr 



h. (41) 



Here the brackets denote the canonical average and trun is the maximum time in 
the run. The result is in Fig. |^. The errors displayed here (and elswhere) cor- 
respond to the canonical (Langevin) average, the microcanonical averages where 
treated as error free. 

In computing the rate we took care that the diffusion {Q'^{t)) is large enough to 
clearly distinguish it from the fluctuations about a classical vacuum (technically 
this is a divergence to be removed by subtraction). These fiuctations can be 
clearly seen in Fig. 0, while in Fig. |^ they correspond roughly to the width ^ 0.5 
of the band. The autocorrelations of this band produce the upward step near 
t = in Fig. H, which can be clearly distinguished from the subsequent linear 



increase with t. In ref. [20| an analytic estimate was made of the contribution of 
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time 

Figure 3: Typical Hamilton evolution of Qhi^) starting with an equilibrium 
configuration at /3 = 13, = 24, v"^ = vl = 0.263. 




time 

Figure 4: The diffusion {Ql{t)) for = 13, = 24, -y^ = -u^ = 0.263, t^un = 
5000. 



13 




200 400 600 BOO 

time 



Figure 5: Same as Fig. ^ after microcanonical averaging with tf, = 4150. 

these fluctuations to the diffusion in SU(2) gauge theory (i.e. the magnitude of 
the step), which roughly fits our data as well. 

The time tb used for the microcanonical averaging (^l]) was varied from zero, 
where t could run up to trun and the diffusion reached about 3-10 (3.5 in the 
case of Fig. ^), to values tb comparable with trun, where the maximum t had to be 
much smaller and the diffusion reached only values of about 1 (0.6 in the case of in 
Fig. 1^). Since such relatively small diffusions correspond to only one unit of AA^^cs 
one may wonder if the diffusion is not dominated by the vacuum fluctuations. 
However, the rates with and without microcanonical averaging were consistent, 
albeit with large errors in the latter case, and furthermore, the contribution of 
vacuum fluctuations can be clearly distinguished, as described above ('the step' 
in Fig. We used /5 = 13 as a test case in the broken phase and similar tests 
were made in the symmetric phase were the rate is larger. Our conclusion is 
that a maximum diffusion of 3 - 10 without microcanonical average is sufficient 
for extracting the rate, using microcanonical averaging with maximum diffusion 
of order 1, as in Fig. |^. This gives us confidence for data points with smaller 
statistics for which microcanonical averaging is necessary for a meaningful result. 

Using the critical values of f3 and v'^ as a starting point we measured the dif- 
fusion {Q'^it)) at various temperatures, volumes and lattice spacings. We mostly 
kept v"^ fixed at v"^ while varying (3, but we also carried out simulations (for (3 = 13 
and 14) in which the deviation of v"^ away from v"^ was calculated according to 
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Figure 6: Results for k = T/L^{awTy as a function of P/Pc- The line connects 
the N = 24: data. 

(PSD, for comparison. The results are recorded in Table |^. Fig. ^ summarizes the 
results for k = T/L^{awT)'^, using Pc to set the scale for the inverse tempera- 
ture. Coming from the lower temperature region we see a jump in the rate at the 
critical temperature, after which k is roughly constant in the high temperature 
phase, hi ^ 1. We should be prepared for the possibility that the flatness of k in 
the high temperature phase may be partly due to the neglect of the renormal- 
ization factor Hg, which is expected to increase rapidly with decreasing (3. The 
falling of the rate at lower /3 < 6 is presumably mainly due to the neglect of this 
finite renormalization Kg. However, there will also be other lattice artefacts in 
this region due to aT getting large (aT = 9/(3). 

The data in Fig. ^ contains tests for volume dependence {N = 8, 20, 24, 
Pc = 12). We see that the > 20 results are volume independent within 
errors, whereas the data for N = 8 show clear finite size effects. The test for 
the dependence on the lattice spacing is less clear cut. We collected data with 
t;2 = corresponding to {Pc,N) = (12,20) and (20,32), with p/pc = 6/12 and 
13/12 (i.e. p = 6,13 for p^ = 12 and p = 10, 21.7 for 4 = 20). The ratio 12/20 
of Pc values corresponds to a ratio (12/20)^ ~ 0.13 in the lattice rate per unit vol- 
ume, aT/N^ = a^r/L^, which is indeed the behavior of the data, approximately. 
Plotting the dimensionless k at the same P/Pc values we see lattice spacing inde- 
pendence within the errors. However, this could be misleading because the left 
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(3 


" c. 


N 


K 




Tc/T 


configs 


4 


12 


24 


0.81(0.05) 


4.19(0.06) 




68 


6 


12 


24 


1.05(0.04) 


2.82(0.04) 




88 


8 


12 


24 


1.09(0.08) 


1.71(0.07) 




71 


10 


12 


24 


1.10(0.06) 


0.85(0.08) 


0.62 


80 


11 


12 


24 


1.09(0.07) 


0.48(0.05) 


0.81 


47 


11.5 


12 


24 


0.92(0.05) 


0.09(0.05) 


0.91 


36 


11.8 


12 


24 


0.66(0.06) 


-0.35(0.09) 


0.96 


27 


12 


12 


24 


0.204(0.009) 


-1.58(0.04) 


1 


44 


12.5 


12 


24 


0.153(0.006) 


-2.04(0.04) 


1.09 


32 


13 


12 


24 


0.141(0.016) 


-2.28(0.11) 


1.18 


33 


14 


12 


24 


0.127(0.009) 


-2.67(0.07) 


1.36 


7 


4 


12 


8 


0.79(0.05) 


4.15(0.06) 




327 


6 


12 


8 


0.77(0.03) 


2.51(0.06) 


_ 


327 


8 


12 


8 


0.276(0.018) 


0.41(0.13) 


- 


93 


6 


12 


20 


1.01(0.03) 


2.78(0.03) 




39 


13 


12 


20 


0.145(0.011) 


-2.25(0.08) 


1.18 


28 


10 


20 


32 


1.09(0.09) 


2.86(0.08) 




30 


21.7 


20 


32 


0.171(0.017) 


-2.08(0.10) 


1.29 


13 


11.5 


12 


32 


0.98(0.06) 


0.14(0.07) 


0.91 


10 


13 


12 


20 


0.163(0.018) 


-2.13(0.11) 


1.08 


17 


14 


12 


24 


0.126(0.010) 


-2.69(0.08) 


1.17 


7 



Table 2: Results for the rate. Here xP' = v'^{Pc, N) as in Table |I], except for the 
two bottom lines, where is related to v1{(3c,N) according to eq. (pSD , keeping 
amn fixed. In the fifth column awTc = l/avr/^c. 

out Kp renormalization factor is different for the two lattice spacings. 

There is another aspect that is relevant here. The lattice spacing dependence 
test should be done for a fixed physical situation, e.g. fixed Tc/T. But using 
P/Pc for the temperature scale is not the same as using T^/T, for the data taken 
at constant v = Vc, since this assumes that the lattice spacing a is constant. 
Changing /5 at constant v implies changing a by eq. (^) (neglecting the running 
of iTiH with T). We can try to calculate Tc/T by calculating ttih/T for each p — v 
pair from (^) and multiplying this by Tc/itlh = 2.15 (from Table |lD. We have 
indicated these Tc/T values in Table ^ for the larger volumes (the perturbative 
formula (|38|) assumes infinite volume). The smaller (3 lead to nonsensical (even 
imaginary) Tc/T, since aT gets too large (recall that aT > 0.75 for f3 < 12). 

Luckily, the constancy of k, in this high temperature region renders the precise 
shift in Tc/T relative to (3 /(3c irrelevant. For the low temperature point the shift 
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0.5 




Figure 7: Plot of k as a function of T^/T. 



in the calculated values of Tc/T = 1.18 {f3 = 13) and 1.29 {f3 = 21.7) away from 
P/Pc = 13/12 ^ 1.08 suggest a slightly larger a-dependence. The two data points 
for which v was calculated using ([38| ) (cf . the two bottom lines in Table ^ , such 
that Tc/T = (3 /(3c, are also consistent within errors. However, the statistics is 
not sufficient to draw a precise conclusion. Plotting the data as a function of 
Tc/T reveals a k, which is remarkably constant also in the low temperature phase 
(Fig. 1^), rather unlike the sphaleron rate. We come back to this in the next 
section. 



6 Conclusion 

We have used dimensional reduction results to pin down the parameters in the ef- 
fective hamiltonian for the classical real time simulations. The difference between 
the effective classical partition function and a more sophisticated dimensional re- 
duction version appears to be numerically insignificant, at least for the critical 
temperature. We conclude therefore that the critical temperature is roughly 
the same as in the DR work in ref. ||2^ (where tuh ~ mw as in our case), 
Tc ~ 2.15 rriH- We tested the rate for lattice spacing dependence in the broken 
as well as in the symmetric phase and found it to be roughly independent of 
a. We take this as evidence that the important parameters in the effective real 
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Figure 8: Comparison of our results for / = ln[r/\^(avy Tc)"^] versus 13/ [3c in the 
broken phase with the sphaleron rate 



time hamiltonian can indeed be deduced from the DR action. We have not taken 
into account yet the finite renormahzation factor /tg (cf. (|40|) ), which may give 
a substantial correction. The neglect of the 'renormahzation of the velocity of 
light' z = ze ^ z-,^ cf. ( |T4| , p8|) is probably insignificant. The errors in our data are 
presumably underestimated, but since we are still in an explorative stage we feel 
this is acceptable. 

In the symmetric phase T > the ratio k = T /V{awT)^ ~ 1 independent 
of T, close to the value 1.1 found in ref. ||2^ for the pure SU(2) gauge theory. 
The rate drops only by a factor of about five when the temperature falls below 
the phase transition. In fact, the numerically obtained rate is much larger than 
the analytic expression suggests: the maximal value of the analytic rate 

is about a factor 650 smaller than our numerical result k ^ 0.20, for T just 
below Tc (cf. appendix B). For comparison we have shown the numerical and 
analytical rate in Fig. |[ The slopes of the numerical rate in the logarithmic 
plot is also much smaller, indicating a rather small effective sphaleron mass. The 
numerical data still seem to be influenced strongly by the phase transition. We do 
not expect, of course, the analytic rate to be reliable near the phase transition, 
but the discrepancy is surprisingly large. Note that taking into account the 
renormahzation factor Kg would enhance the difference. 

Another puzzling indication for the discrepancy is the following. If the rate 
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in the low temperature phase were given by the analytic expression, we would 
most probably not have seen a single transition in the time span in figure 0, for 
the following reasoning. If the transitions with ANq ^ 1 are as clear as in Fig. ^ 
the rate may be calculated as F = n/t, where n is the number of transitions that 
have occured in time t, provided n,t ^ oo. Using n = 3 and t = 1.4 * 10^ gives 
the rough estimate ln[F/L'^(avi/Tc)^] = —4.8, with an perhaps error of a factor of 
2 because n is small. This estimate corresponds to a point in figure ^ near the 
linear extrapolation of the numerical data to /5//5c = 3, far above the analytical 
curve there. The temperature in Fig. ^is rather low, T/Tc ~ T/2mu ^ 1/3, and 
one may doubt the applicability of dimensional reduction in this case. However, 
the analytic calculation is also based on dimensional reduction and it makes sense 
to make the comparison with the numerical results. 

Since the dependence of the rate on the sphaleron energy Eg is exponen- 
tial, lattice artefacts in this quantity could give large systematic errors. Recent 
calculations of Eg |^ have shown that the difference AE^ = Es{a) — Es{0) 
{a = lattice spacing) is negative, which would enhance the rate on the lat- 
tice. Although we did not find a significant lattice spacing dependence near 
the phase transition, it is of interest to estimate the systematic error. Ref. 
3^ found that AB = awAEs/2mw ~ ca^m^, with c = —0.12. Identify- 



ing this m\Y with the temperature dependent screening mass mw{T), we can 
estimate exp{—AEs/T) ^ exp[\c\a'^rn'^{T){2mw(T)/T){l/aw)]- In the broken 
phase mw(T) decreases roughly from Tc/2 at T = to Tc/3 at T = Tc It 
follows that 2mw{T)/T < 2 in the region Tc/2 < T < T^, and with aT^ = 3/4 for 
(3c = 12, we also have a?myY{T) < a^T^/^ = 9/64. Putting things together this 
gives exp{—AEs/T) < 2.6 in the region Tc/2 < T < which amply includes 
our data. Such a large systematic effect on the rate would be a substantial error 
indeed, but it cannot explain a factor > 650. 

Assuming our results to hold up to future scrutiny, they suggest that many 
non-sphaleron processes are contributing to the rate in the broken phase. A 
classification of such configurations is given in [0. The danger of erasure of 



a surplus of baryon number after the electroweak phase transition in the early 
universe, which is already a possibility with the usual sphaleron rate, is of course 
only magnified by these results. 
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A Canonical formalism for lattice gauge fields 

For completeness we record here some details on the canonical formalism with lat- 
tice gauge fields using continuous time @]. Our main objective is to understand 
the measure YlidUdE) used in the classical partition function in relation 
to the canonical measure ]^((iy4(in). As an example we consider a one-link la- 
grangian of the form 

L = ^ Tr f/f/t = -1^ (42) 

where L'" is a group element with coordinates and Gap = 2Tr [{dU / dA^) {dUydA^] 
is a metric on group space. The canonical conjugate to A" is given by 



Iia = ^GapA^, (43) 

and the hamiltonian is 



H = ]^zg^G'^PTlJlp, (44) 

where G"^ is the inverse of Gap- G^pG^"' = 51. In the quantization of the system 
it is convenient to introduce left and right handed vielbeins according to 

dU 

Va{R) = -zt/t__ = y„^(i?)t^, (45) 

= iU^ = Vap{L)tp, (46) 

where the ta are generators in the fundamental representation of the gauge group. 
We assume these to be normalized as It tat p = 6013/2 (for SU(2) ta = Ta/2). In 
terms of the Va/3 and their inverse V^^ (i.e. VapV^i^ = 5^ for R and L, respectively), 
one introduces the £"s as 

Ea{R) = Vi{R)Ufs, (47) 

and similar for L. The Ea{L) and Ea{R) are the generators of left and right 
translations in the sense that 

[Ea{L),U]=taU, [Ea{R),U] = Uta, (48) 

follow from the canonical commutation relations [A", 11^] = Since G"^ = 
V"{L)V'^{L) = V'^{R)V/^{R), the hamiltonian can be rewritten as 

H = ^ zg^Ea{L)Ea{L) = ^ zg^Ea{R)Ea{R), (49) 
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which provides an attractive prescription for the operator ordering problem in 
the quantum theory. 

Returning to the classical theory, the canonical partition function Z^. is given 
by the integral of eyi]i{—H /T) over phase space with volume element ^^((iA^dnQ,). 
As usual, this volume element is invariant under transformations of coordinates 
y4°. It is, however, generally more transparant to trade the II's for the i?'s. 
Then the invariant volume element on group space dU arises naturally from 

= V^p{L)Ep{L), 

W^dA'^dli^) = I det V{L)\ YlldA^'dEo.iL)] = Vd^ Y[[dA"dE^{L)] 

a a a 

= const. dU YldEc,{L), (50) 

a 

and similar for R. 



The conventions used in ref. []T8[ can be reached from those in sect. 3 by 
first absorbing z in the time scale and then writing L = {A/g'^)L. This L = 
{1/A)TtUW is equal to L with z = A/g"^ and going through the canonical formal- 
ism again gives H = 2E'^ in terms of the new canonical -E's. The E^ are related 
to these by E^ = 2Ea. 

B Analytic rate 

We give here some details on the analytic form (0,^ for the rate. The function 
/(-^3/5'i) in (D is given by 

where p = \uj^\/2mw with oj^ the unstable sphaleron eigenvalue, B = awEs/2mw, 
TVtr and A^rot are zero mode factors and S is the ratio of fluctuation determinants. 
We treat the temperature dependence of various quantities in a customary ap- 
proximation, in which p, B, A/tr, AArot and S are independent of T and depend 
only on Xs/g^ ~ 8m|^/m^. All temperature dependence is then in mw which is 
approximated as 

mw (T) =mw{0) y^l - T^/T^ (52) 

For the case rriH = mw we have p ~ 0.73, B ^ 1.8, A/tr ~ 7, A/'rot ~ 12 [|] and 
6 ~ exp(— 9.64) 0. The temperature dependence in eq. (^ now enters through 
the factor a;^exp(— x), x = 2Bmw{T)/awT, which has a maximum of ~ 751 at 
X = 7, or T/Tc ~ 0.99. The rate at this maximum is given by k ~ 3.1 x 10""^, 
which is a factor of about 650 smaller than the value ~ 0.20 found numerically 
for T just below the phase transition. 
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